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ABSTRACT 

The prolific magnetar SGR 1900+14 showed two outbursts in the last decade and has been closely 
monitored in the X-rays to track the changes in its radiative properties. We use archival Chandra and 
XMM-Newton observations of SGR 1900+14 to construct a history of its spectrum and persistent X- 
ray flux spanning a period of about seven years. We show that the decline of its X-ray flux in these two 
outburst episodes follows the same trend. The flux begins to decline promptly and rapidly subsequent 
to the flares, then decreases gradually for about 600 days, at which point it resumes a more rapid 
decline. Utilizing the high quality spectral data in each epoch, we also study the spectral coevolution 
of the source with its flux. We find that neither the magnetic field strength nor the magnetospheric 
properties change over the period spanned by the observations, while the surface temperature as well 
as the inferred emitting area both decline with time following both outbursts. We also show that 
the source reached the same minimum flux level in its decline from these two subsequent outbursts, 
suggesting that this flux level may be its steady quiescent flux. 
Subject headings: pulsars: individual (SGR 1900+14) — X-rays: bursts 



1. INTRODUCTION 

Soft Gamma Repeaters (SGRs) and Anomalous X-ray 
Pulsars (AXPs) belong to a class of objects called magne- 
tars - neutron stars whose X-ray emission is likely to be 
powered by the decay of their extremely strong magnetic 
fields (Duncan & Thompson 1992; Thompson & Duncan 
1996; Thompson, Lyutikov & Kulkarni 2002). All seven 
confirmed SGRs and six out of seven confirmed AXPsQ 
have emitted energetic bursts of X-rays/soft gamma rays. 
Burst active episodes of magnetars last anywhere from 
few hours to months. During their bursting activity, 
magnetars also exhibit remarkable temporal and spec- 
tral changes in their persistent X-ray output. A detailed 
description of SGRs and AXPs can be found in Woods 
& Thompson (2006) and Mereghetti (2008). 

In the last seven years, new magnetar candidates have 
emerged, most prominently through transient outbursts 
(e.g., XTE J1810-197, Halpern & Gotthelf 2005; CXO 
J164710.2-455216, Israel et al. 2007; SGR J1833-0832, 
G6gii§ et al. 2010). These sources were too dim to be 
detected in X-rays during their quiescent phases (usually 
below our detection sensitivity), but their X-ray fluxes 
increased by up to few hundred times as they entered 
their outburst episodes. In addition, known magnetar 
sources also exhibit variations (triggered by bursting ac- 
tivity) in their persistent flux, although typically not as 
dramatic as those in the transient systems (e.g.. Woods 
et al. 2004). These high X-ray luminosities were in- 
strumental in probing the outburst mechanism of mag- 
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netars (Ozel & Giiver 2007; Giiver et al. 2007; Ng et 
al. 2010). In contrast, the low but persistent flux level 
of most magnetars requires long term (5-10 years) mon- 
itoring to understand (burst-induced) changes in their 
emission properties (Dib, Kaspi, & Gavriil 2009). 

SGR 1900+14 has been one of the most prolific SGRs: 
it was discovered in 1979 (Mazets, Golenetskij & Guryan 
1979), and was detected in a bursting mode again in 1992 
(Kouveliotou et al. 1993). In May 1998 the source en- 
tered a major outburst episode that lasted about eight 
months and included the giant flare on 1998 August 27 
(Hurley et al. 1999). ASCA and RXTE observations 
prior to and during the 1998 activation led to the dis- 
covery of its 5.16 s spin period (Hurley et al. 1999), 
its spin-down rate of ~10~^^ s/s and magnetic field of 
2 - 8 X 10^4 G (Kouveliotou et al. 1999), and thus to 
the confirmation of its magnetar nature. The source re- 
sumed a high level of activity in April 2001 (Guidorzi et 
al. 2001; Kouveliotou et al. 2001) and again in March 
2006 (Vetere et al. 2006). As its burst active phases are 
well separated from each other, SGR 1900+14 is an excel- 
lent source to investigate radiative changes both during 
bursting behavior as well as in burst quiescence. 

The first major enhancement in the persistent X-ray 
flux of SGR 1900+14 was observed at the onset of the 
1998 August 27 giant flare: the flux increased by a fac- 
tor of ^700 with respect to its level before the activa- 
tion (Woods et al. 2001). A detailed spectral analysis 
of the (much longer) flux decay period was not possible 
due to the lack of continuous monitoring observations 
with imaging instruments during the decay phase. Es- 
posito et al. (2007) analyzed nine pointed BeppoSAX 
observations of SGR 1900+14 spanning about five years 
from May 1997 to Aprfl 2002 and noted that the (2-10 
keV) fiux measured in the last pointing faded signifi- 
cantly with respect to earlier observations. Following the 
April 2001 activation, the Chandra X-ray Observatory 
and XMM-Newton observed SGR 1900+14 at numerous 
occasions, establishing a valuable dataset for understand- 
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ing the long term behavior of this source in, particular, 
and of magnetars, in general. 

In this paper, we make use of all archival Chan- 
dra and XMM-Newton observations to construct the 
persistent X-ray flux temporal and spectral history of 
SGR 1900+14 spanning about seven years following the 
April 2001 activation. In the next section we introduce 
the Chandra and XMM-Newton observations used in this 
study. In Section 3, we present the results of the spec- 
tral analysis and show that the magnetic field strength 
and the magnetospheric properties remain stable over 
the period spanned by the observations, while the sur- 
face temperature and the inferred emitting area both de- 
cline with time following both outbursts. We also show 
that the source flux shows the same trend in its decline 
from outburst in both episodes and ultimately reaches 
the same minimum flux in both cases. We discuss the 
implications of these results in Section 4. 

2. OBSERVATIONS AND DATA ANALYSIS 

Between 2001 April 22 and 2008 April 8, SGR 1900+14 
was observed 13 times with Chandra and XMM- 
NewtonTahle [1] lists the log of these pointed X-ray ob- 
servations. We describe below the details of our data 
reduction. 
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TABLE 1 
Log of SGR 1900+14 Observations 

Observatory Observation ID Observation Exposure Time 







Date 


(ks) 


Chandra 


2458 


2001 Apr 22 


20.1 


Chandra 


2459 


2001 Apr 30 


18.8 


Chandra 


3858 


2002 Nov 6 


48.0 


Chandra 


3862 


2003 Feb 18 


25.1 


Chandra 


3863 


2003 Jun 2 


25.6 


Chandra 


3864 


2003 Oct 18 


25.3 


XMM-Newton 


0305580101 


2005 Sep 20 


20.2 


XMM-Newton 


0305580201 


2005 Sep 22 


18.7 


Chandra 


6709 


2006 Mar 29 


40.0 


XMM-Newton 


0410580101 


2006 Apr 1 


13.4 


Chandra 


7593 


2007 Jun 24 


12.2 


Chandra 


8215 


2007 Nov 21 


12.6 


XMM-Newton 


0506430101 


2008 Apr 8 


18.7 



2.1. Chandra 

There are a total of nine Chandra observations per- 
formed between 2001 Aprfl 22 and 2007 November 21. 
These were all performed using the Advanced CCD Imag- 
ing Spectrometer (ACIS) in continuous clocking (CC) 
mode. We selected rectangular source regions centered 
at the position of SGR 1900+14 with dimensions 8" x 2". 
Our background regions were selected with similar rect- 
angular sizes from source-free regions on the collapsed 
CC mode image. We calibrated the Chandra observa- 
tions using the CIAO version 4.2 and the CALDB version 
4.2. 

We extracted the X-ray spectra following the standard 
Chandra data analysis threads with the psextract tool. 
We then used the mkacisrmf and mkarf tools to create 
the detector and ancillary response files, respectively. Fi- 
nally, we re-binned the spectra such that each energy bin 



would contain at least 50 counts, to decrease the formal 
errors. 

2.2. XMM-Newton 

We analyzed the X-ray data obtained with the EPIC- 
pn detector between 2005 September 20 and 2008 April 
8. In all observations, the EPIC-pn was operating in full 
frame mode. The calibration of the data was performed 
with the Science Analysis Software (SAS) version 9.0.0 
and the latest available calibration files as of February 
2010, using the task epproc. 

We extracted the X-ray source spectra by accumulat- 
ing events from a circular region with a radius of 32" 
centered at SGR 1900+14. The background regions were 
selected on the same chip from source free regions with 
a typical radius of 50". We then used the rmfgen and 
arfgen tools of SAS to generate detector and ancillary 
response files and re-binned the X-ray spectra such that 
each energy bin contained at least 50 counts. 

3. SPECTRAL ANALYSIS AND RESULTS 

We used XSPEC vl2.5.1n (Arnaud 1996) to analyze 
all spectra. We assumed a gravitational redshift correc- 
tion of 0.306, corresponding to a neutron star mass of 
1.4 Mq and a radius of 10 km. We performed the fits 
in the 0.8—6.5 keV range, where the lower energy bound 
is set by the source flux and the higher energy limit is 
determined by the non-thermal hard X-ray component 
(Gotz et al. 2006) affecting the soft X-ray spectra. 

We fit the X-ray spectra using the Surface Ther- 
mal Emission and Magnetospheric Scattering (STEMS; 
Oliver, Ozel and Lyutikov 2006) model. STEMS is based 
on the radiative equilibrium atmosphere calculations pre- 
sented in Ozel (2001, 2003) but also includes the effects of 
magnetospheric scattering on the photons emitted from 
the neutron star surface as calculated by Lyutikov & 
Gavriil (2006). The model parameters consist of the sur- 
face magnetic field strength, B, surface temperature, T, 
the magnetospheric scattering optical depth, r, and the 
velocity of the particles in the magnetosphere, /3 = v/c, 
where c is the speed of light. 

We first fit all 13 spectra simultaneously using STEMS 
and taking into account the effect of interstellar absorp- 
tion. Wc allowed the spectral parameters of all spectra 
to vary individually. While we did not specify its value 
a priori, we forced the hydrogen column density to be 
the same between observations. We found that the mag- 
netic field strength, the scattering optical depth, and the 
average particle velocity did not vary significantly be- 
tween different observations. We measured a hydrogen 
column density of (2.36 ± 0.05) x 10^^ cm~^, assuming 
solar abundances (Anders & Grevesse 1989). We, there- 
fore, performed all further STEMS fits fixing the hydro- 
gen column density at the above value and forcing the 
parameters B, t, and /3 to be constant between different 
observations. We still allowed the temperature and the 
model normalization to vary individually for all spectra. 
We obtained good fits to all 13 spectra, with a x^/dof 
= 2370/2438. In this combined fit, we found a magnetic 
field strength of B = (5.0 ± 0.3) x lO^^ G, a scattering 
optical depth of r = 8.7 ± 0.9 and a particle velocity of 
(3 = 0.37 ± 0.01. Below, we discuss the time evolution 
of the temperature T and the source flux based on these 
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Fig. 1. — Upper panel: two X-ray spectra of SGR 1900+14 ob- 
served with the XMM-NewtonEPlC-pn (upper data set in black) 
and the Chandra/ ACTS in CC mode (lower data set in red). Solid 
lines are the best fitting STEMS model curves. Lower panel; resid- 
uals of the spectral fits. 

combined fits. Note that errors reported throughout the 
paper are 1 a. 

We also investigated the potential effects of the cross 
calibration between the Chandra/ ACIS in CC mode 
and the XMM-Newton/EPlC-pn on the determination 
of spectral parameters. To this end, we used two data 
sets, observation ID 7593 with Chandra and 0410580101 
with XMM-Newton which were selected because they oc- 
curred at comparable source fluxes. We fit both spectra 
using STEMS as well as the empirical blackbody plus 
power law model, accounting for interstellar absorption 
in both cases. In this analysis, we again forced the ab- 
sorption parameter to remain constant for both spectra 
as it is not observed to vary over time. In Figure [U 
we present the two spectra along with the best fitting 
STEMS models. We found all parameters of both con- 
tinuum models to be consistent between the Chandra and 
XMM spectra to within 1 a errors. These results ensure 
that the use of two different instruments does not intro- 
duce any systematic biases in the joint spectral analysis. 

We present in Figure [2] the history of the unabsorbed 
X-ray fiux of SGR 1900-1-14 in the 0.8-6.5 keV band. 
The first Chandra observation took place four days after 
the intermediate event (Feroci et al. 2002). The source 
flux declined rapidly in the period soon after this event, 
dropping from (1.34 ± 0.02) x 10~^^ erg cm^^ s^^ to 
(1.16±0.02) X 10"-'^-'^ erg cm~^ s~^ in only about six days. 
An even faster decline was seen in the contemporaneous 
iJeppoS'^X observations: the source fiux in the 2—10 keV 
band was (3.1 ±0.3) x 10~" erg cm^^ g-i 2OOI April 
18 and declined to (1.06 ± 0.04) x 10"" erg cm^^ 
in 11 days (Esposito et al. 2007). We observed a similar 
trend following the 2006 reactivation of the source as 
its fiux drops from (8.6 ± 0.1) x 10~^^ erg cm~^ s~^ to 
(7.9 ± 0.1) X 10^^^ erg cm^^ s^^ in three days. 

We further studied these long term persistent fiux vari- 
ations of SGR 1900-1-14 as follows: we determined the 
relative times of the first eight observations (which took 
place before 2006) with respect to the onset of the April 
2001 outburst and those of the remaining five observa- 
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Fig. 2.— Unabsorbed flux history of SGR 1900-1-14 in the 0.8-6.5 
keV range. The two arrows indicate the onset of the April 2001 and 
March 2006 reactivations of the source, respectively. The calendar 
dates on top of the figure are in the YY/MM/DD format. 

tions with respect to the onset of the March 2006 ac- 
tivation. We present in Figure [3] the unabsorbed flux 
as a function of relative time since each respective out- 
burst onset. We found that in both outbursts the source 
flux declined rather gradually until about 600 days after 
the onset and exhibited a sharper decline trend beyond 
- 600 days. We also found that SGR 1900+14 was at 
its lowest X-ray fiux level of 6.8 x 10~^^ erg cm~^ s~^ 
before the March 2006 reactivation. Following the source 
re-brightening after the 2006 outburst, the flux reached 
that level again in the last pointed observation. 

To better understand the nature of these flux variations 
in SGR 1900+14, we investigated a possible correlation 
between the flux and the only varying STEMS parame- 
ter, i.e., the surface temperature. We find that the flux 
and surface temperature are indeed correlated, with a 
Spearman's rank order correlation coefficient r = 0.72. 
The probability of obtaining such a correlation with a 
random data set is P = 0.005. In Figure IH we show the 
history of the surface temperature of the neutron star 
as well as its long term fiux behavior. The decline in 
the surface temperature alone (assuming a single tem- 
perature) cannot account for the observed flux decline. 
In particular, the average surface temperature dropped 
from ~ 0.56 keV to 0.53 keV as measured on 2001 
April 22 (MJD 52021) and 2005 September 22 (MJD 
53635), respectively, which corresponds to a flux decline 
of about 25% (if the emitting surface area remains con- 
stant). However, the source flux declined by about 97% 
between these two epochs, mostly soon after the rise; the 
remaining decline over the last 1100 days of the observ- 
ing period was only about 33%, while the temperature 
decreased by about 6%. 

We also investigated the long term behavior of the sur- 
face emitting area inferred using the STEMS normal- 
ization and surface temperature. We found that, for a 
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Fig. 3.— Unabsorbed flux of SGR 1900+14 in the 0.8 - 6.5 keV 
range vs. relative time since the onset of the two bursting episodes: 
2001 April (triangles) and 2006 March (squares). 
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Fig. 5.— Pulse profile of SGR 1900+14 during the Chandra ob- 
servations on 2001 April 22 (MJD 52021). The pulse phase intervals 
used as pulse peak and pulse minimum spectral accumulations are 
shown with hatched and dotted regions, respectively. 
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Fig. 4. — Evolution of the neutron star surface temperatures 
(squares) obtained by fitting each spectrum with STEMS, and the 
corresponding unabsorbed flux (diamonds). The vertical dashed 
lines indicate the onsets of the 2001 and 2006 outbursts, respec- 
tively. Arrows at the lower end of the flgure indicate observations 
used in further phase resolved spectral analysis (see the text). Cal- 
endar dates on top are in the YY/MM/DD format. 



uniform surface temperature, the radius of the emitting 
region is 4.8±0.2 km soon after the 2001 outburst on- 
set (assuming a distance to SGR 1900+14 of 13.5 kpc; 
Vrba et al. 2000). It then dechned to 3.9±0.2 km and 
remained fairly constant until 2005. The radius emit- 
ting surface went up again to 4.4±0.1 km following the 
2006 outburst but quickly fell to a marginally higher con- 
stant level of 4.1±0.2 km. This change, however, could 
be the result of neglecting temperature inhomogeneities. 
Note also that the errors in the surface temperature and 
normalization, and, consequently, in the radius of the 
emitting area, are correlated. 

Finally, we performed a coarse phase resolved spectral 
analysis using five Chandra observations (indicated by 
arrows in Figure |4]) to check whether there are signifi- 
cant surface temperature variations over the spin phase, 
which may affect or bias the spectral determination of 
the phase-averaged temperature. For each of the selected 
pointings, we obtained a pulse peak spectrum (spanning 
0.30 of the spin phase during the pulse peak) and a pulse 
minimum spectrum (spanning 0.30 of the spin phase dur- 
ing the pulse minimum) using accurate and contempora- 
neous pulse period reported in Mcrcghetti et al. (2006). 
We show in Figure [5] the pulse profile and pulse phase in- 
tervals within which the pulse peak and pulse minimum 
spectra were obtained in the observation on 2001 April 
22 (MJD 52021). We fit the X-ray spectra in the 0.8-6.5 
keV using STEMS with surface magnetic field strength 
and magnetospheric parameters fixed at their phase av- 
eraged values. In Table [H we list the temperatures and 
model normalizations (i.e., an indicator of intensity) for 
the five selected observations. We find that the surface 
temperature remains constant over the spin phase of the 
source in the majority of these observations. 
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TABLE 2 

Spectral Fit Results of the Phase-resolved Analysis 



Obs Date 


Pulse Peak 


Pulse Minimum 


Peak Normalization 


Min Normalization 


(MJD) 


Temperature (keV) 


Temperature (keV) 


(10-"') 


(10-"') 


52021 


0.57±0.01 


0.54±0.01 


2.2±0.3 


1.8±0.1 


52584 


0.56±0.01 


0.55±0.01 


1.6±0.1 


1.2±0.1 


52930 


0.53±0.01 


0.52±0.01 


0.57±0.05 


0.48±0.05 


53823 


0.53±0.01 


0.52±0.01 


0.26±0.03 


0.19±0.02 


54275 


0.55±0.01 


0.47±0.01 


0.37±0.03 


0.5±0.2 



4. DISCUSSION 

We found that SGR 1900+14 exhibits a monotonic but 
non-steady flux dechne following the X-ray brightening 
during its 2001 and 2006 outbursts. In both outbursts, 
its flux dropped rapidly within a lew weeks after the 
onset of the outburst and then at a slower rate for ap- 
proximately 600 days. After this period, the flux declines 
again at a much faster rate. A similar decay trend was 
also seen in SGR 1627-41 (Kouveliotou et al. 2003). The 
fact that we observe the same trend in successive out- 
bursts from the same source suggests that the long-term 
effects of outbursts are not stochastic but reproducible. 

It is clear that magnetar bursting activity leads to long- 
term flux enhancements and that the additional energy 
powering these enhancements is stored in the crust as 
heat. This heat comes most likely from the energy re- 
leased in the crust during the bursts, or from the energy 
deposited in the crust by the bombardment with mag- 
netospheric particles during the burst. The crust, then, 
reradiates this additional heat over a timescale of a few 
years. 

Despite the correlation between the flux and the sur- 
face temperature, the variation in the surface tempera- 
ture alone does not account for the decline in the flux, 
but it also requires a change in the emitting area over 
time. This can perhaps be understood if the crust is 
heated inhomogeneously, as would be expected if the ini- 
tial heating episode is due to magnetic energy release. 
Moreover, because the thermal resistance of the crust 
is dominated by the uppermost layers, where the heat 
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conductivity is strongly affected by the magnetic fleld, 
heat coming from the deeper layers of the crust could 
reach the surface unevenly. This leads to both uneven 
heating and uneven cooling, which may affect the total 
inferred emitting area. Furthermore, over-time, because 
the observations are carried out in a limited energy range, 
cooler parts of the crust may fall out of the observed en- 
ergy band faster than the hotter regions, again reducing 
the inferred emitting area. Thus, the observed source 
flux would decline both due to a decline in temperature 
as oc T'^ and due to a decline in the emitting area. 

We flnd on three different occasions that the X-ray flux 
of SGR 1900-H14 was as low as 6.8 x lO'^^ gj-g cm'^ s-^, 
suggesting that this level may correspond to the source 
persistent X-ray flux in the absence of burst induced en- 
hancements. The current detection thresholds of imaging 
instruments are ~10~^'^ erg cm"^ s~^; it is, therefore, 
possible that the persistent flux levels of the so-called 
transient magnetars are much lower than those of the 
always detectable magnetars. If that is indeed the case, 
the magnetar engine that is responsible for a source per- 
sistent quiescent X-ray emission would seem to power a 
broad flux range of < 10^^'^ to ~ 10^^" erg cm-^ s-^. 
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